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Abstract 

We analyze the nearest neighbor spacing distributions of low-lying 2 + levels of 
even-even nuclei. We grouped the nuclei into classes defined by the quadrupole 
deformation parameter (/? 2 )- We calculate the nearest neighbor spacing distributions 
for each class. Then, we determine the chaoticity parameter for each class with the 
help of the Bayesian inference method. We compare these distributions to a formula 
that describes the transition to chaos by varying a tuning parameter. This parameter 
appears to depend in a non-trivial way on the nuclear deformation, and takes small 
values indicating regularity in strongly deformed nuclei and especially in those having 
an oblate deformation. 

1. Introduction 

There are two main ways to investigate nuclear structure; one of them is the 
detailed study of individual nuclei. This method tests different nuclear models, and 
gives us the internal structure of the nucleus under investigation. However, this 
method does not allow constructing a correlation relation between different nuclei. In 
addition, we are in need to investigate a large number of studies to understand nearly 
2500 discovered nuclei. The second way is the cumulative study, in which one 
attempts to classify nuclei in terms of a small number of parameters that can be 
related to nuclear structure. This method finds the correlations between different 
nuclei and helps us to understand the nuclear structure evolution through the nuclear 
chart. 

Several papers have been published [see, e.g. 1,2,3,4,5], whose purpose is to find 
the extent of chaos in nuclear dynamics. Random matrix theory (RMT) has gained a 
great success in the description of the fluctuation properties of spectra of quantum 
systems. It models a chaotic system by an ensemble of random matrices subject only 
to symmetry constraints. Systems conserving time reversal, such as the atomic 
nucleus are described by the Gaussian orthogonal ensemble (GOE) [6]. 

The nearest neighbor spacing (NNS) distribution occupies an important place in 
the statistical theory of spectra. Unfortunately, RMT does not provide a closed form 
expression for the NNS distribution, which is suitable for the analysis of experimental 
data. However, Wigner proposed an approximate expression, which is exact in the 
case of 2 x 2 matrices. The NNS distribution of nuclear spectra at neutron threshold 
energies is in good agreement with the Wigner distribution [7]. The study of level 
statistics at low-lying excitation energies requires complete (few or no missing levels) 
and pure (few or no unknown spin-parities) level schemes. Unluckily complete and 
pure level schemes are only available for a limited number of nuclides. Therefore, to 
improve the statistical significance of the study, spacing distributions for several 
sequences of levels (each sequence has the same spin parity) from different nuclides 
have been combined. 

Such statistical studies concluded that the NNS distribution of nuclei in the ground 
state region is intermediate between the Poisson distribution expected for regular 
systems and the Wigner distribution for chaotic ones. In addition, light nuclides 
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showed a behavior close to chaotic near the ground state, while heavier nuclides 
seemed to be more regular. This can be understood from the fact that the fluctuation 
properties of the spectra given by the nuclear shell model (which well describes light 
nuclei) agree with the GOE. 

In the present study, we investigate the effect of nuclear deformation -as a measure 
of nuclear structure- on the chaoticity of even-even nuclei. Due to serious statistical 
limitation of the experimental data basis, earlier results yielded only qualitative 
information concerning the effect of deformation. Some efforts have focused on the 
symmetries of the IBA, but again the statistics were too limited to reach definite 
conclusions. In this paper we wish to establish a direct relation between deformation 
and nuclear level statistics by parameterizing the available experimental data via 
quadrupole deformation parameter. 

We focus on the 2 + states because of their abundance in even-even nuclei. The data 
set and the deformation parameter used to classify the nuclei are described in Sec. 2. 
Section 3 describes the technique of the analysis, while the results are given in Sec. 4. 

2. Data Set 

In this section, we describe our choice of levels in even-even nuclei. We recall that 
the collective motion is interpreted as vibrations and rotations of the nuclear surface 
in the geometric collective model first proposed by Bohr and Mottelson [8], where the 
nucleus is modeled as a charged liquid drop. The moving nuclear surface may be 
described by an expansion in spherical harmonics with time-dependent shape 
parameters as coefficient. The quadrupole deformation seems to be the most 
important collective excitations of the nucleus. For axially symmetric nuclei, the 
nuclear radius can be written as, 

R{e,4>)=R m [\+p?n(°M, (2) 

where P2 is the quadrupole deformation parameter. Positive and negative P2 values 
correspond to prolate and oblate shapes respectively. The deformation parameter 
suffers from the difficulty to distinguish between static deformation and the dynamic 
amplitude of the quadrupole vibration in soft spherical nuclei. It can be obtained only 
by a detailed analysis of spectra and transition probabilities. Despite this difficulty, we 
draw a crude picture of the dependence of nuclear level statistics on nuclear 
deformation in a direct way. 

The quadrupole deformation parameters P2 are taken from macroscopic- 
microscopic calculations [9]. In order to obtain a sufficient number of energy levels 
within each interval, we take the absolute value of P2. The Z dependence of 
deformation is not taken into account. We compare results by analyzing the data set in 
terms of the experimental P2 values, as deduced from the B(E2; + — > 2 + ) values [10]. 

The data on low-lying 2 + levels of even-even nuclei are taken from Endt [11] for 
22 < A < 44, and from the Nuclear Data Sheets until March 2005 for heavier nuclei. 
We consider nuclei in which the spin-parity J 1 assignments of at least five consecutive 
levels are definite. In cases where the spin-parity assignments are uncertain and where 
the most probable value appeared in brackets, we accept this value. We terminate the 
sequence in each nucleus when we arrive at a level with unassigned J 1 , or when an 
ambiguous assignment involved a spin-parity among several possibilities, as e.g. J 1 = 
(2 + , 4 + ). We chose one of the suggested assignments, when only one such level 
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occurred in the sequence, and was followed by several definitely assigned levels 
containing at least two levels of the same spin-parity, provided that the ambiguous 
level is found in a similar position in the spectrum of a neighboring nucleus. 
However, this situation has occurred for less than 5% of the levels considered. In this 
way, we obtained 1 132 levels of spin-parity 2 + belonging to 150 nuclei. 



3. Method of Analysis 



For statistical studies using RMT, one should have a spectrum of unit mean level 
spacing. This is obtained by fitting a theoretical expression to the number N(E) of 
levels below excitation energy E; this process is called unfolding. The expression 
used here is the constant-temperature formula (4), 

N{E) = N + aqp(^A. (4) 

V 1 ) 

The three parameters No, Eo and T obtained for each nucleus vary considerably with 
mass number. Nevertheless, all three show a clear tendency to decrease with 
increasing mass number. A detailed account of the method of analysis used in the 
present work has been given in Refs. [12, 13], and references therein. 

The nuclear states are characterized by their invariance under time reversal and 
space rotation, which can be represented by the Gaussian orthogonal ensemble (GOE) 
of random matrices. The NNS distribution of levels of the GOE is well approximated 
by Wigner's distribution [7] 

* ( ~-A- (5) 
4 J 



Pw (s ) = —s exp 



Here, s is the spacing of neighboring levels in units of the mean level spacing. For 
integrable systems, the NNS distribution is generically given by the Poisson 
distribution, 

P P (s) = exp(-s). (6) 
The key ingredient of the present analysis is the assumption that the deviation of the 
NNS distribution of low-lying nuclear levels from the GOE statistics is caused by the 
neglect of possibly existing conserved quantum numbers other than energy, spin, and 
parity. A given sequence S of levels can then be represented as a superposition of m 
independent sequences Sj each having fractional level density fj, with j = 1, . . .,m, and 

with <fj < 1 and £" =1 // = 1 . We assume that the NNS distribution Pj(s) of Sj obeys 

GOE statistics. The exact NNS distribution P(s) of this superposition has been given 



in Ref. [6]. It depends on the (m-1) parameters fj, j = 1, 



,m-L In [14], this 



expression has been simplified by observing that P(s) is mainly determined by short- 
range level correlations. This reduces the number of parameters to unity and the 
proposed NNS distribution of the spectrum is 



P{s,f) = 



ns 



1-/4-/(0.7 + 0.3/) — 



:exp 



-(!-/>-/ (0.7 + 0.3/) 



ns 



2 \ 



(7) 



which depends on only one parameter, the mean fractional level number 

/ = '^ l "j_ l fj 2 for the superimposed sub-spectra. This quantity will eventually be used 

as a fit parameter. For a large number m of sub-spectra, / is of the order of \lm. In 
this limit, P{s, f) approaches the Poisson distribution P(s, 0) = Pp(s). This expresses 
the well-known fact that the superposition of many GOE level-sequences produces a 
Poissonian sequence. On the other hand, when / — > 1 the spectrum approaches the 
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GOE behavior. Indeed, P(s, 1) coincides with the Wigner distribution (5), this is 
expected as the system then consists of a single GOE sequence. This is why /is called 
the chaoticity parameter. 

In this study, we use the parameter / instead of empirical Brody's interpolation 
formula P(s,&>) = a>(y + \)s r exp(-&>s r+1 ), where y is a fitting parameter and 



r+2 



y+1 



[15]. Since, formula (7) has a theoretical basis because / represents the 



area of phase space of the wave function occupied by the chaotic dynamics. We plot a 
relation in figure 1 between the variation of parameter / and the parameter co, by 
calculating the chi-square fit between the NNS distribution in Eq. (7), and Brody 

5 

distribution. % 2 = \{P( S J )—P(s,co)) 2 ds . By defining the parameter/ we calculate 

o 

the co parameter. We observe nearly a linear relationship between the two parameters, 
that strength our choice of parameter / 
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Fig. 1: A relation between a parameter / and co of Brody distribution. 



A common method to determine the parameter / is via the best fit obtained by chi- 
square. But working with histograms of experimental data, there is always the issue of 
selection of bin size. Different choices of the bin width, particularly when the widths 
vary, can lead to quite different visual conclusions. In order to avoid this difficulty, 
Bayesian inference [16] is used to determine the parameter / The Bayesian analysis is 
independent of the bin size; it deals with the spacings directly [13]. The Bayesian 
analysis yields the best-fit value of the chaoticity parameter/ and its error for each 
NNS distribution. For a given sequence of spacings s = (si, si, ... sn), the joint 
probability distribution P(s\f) of these spacings, conditioned by the parameter /, is 
given by 



P(s\f) = YlP(s i ,f), 



(8) 



with P(si, J) given by Eq. (7). Bayes' theorem provides the posterior distribution 

WW = ™, (9) 
M(s) 

of the parameter / given the events s. Here, ju(f) is the prior distribution and 
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M( S ) = \ 1 o P(s\f) J u(f)df, 



(10) 



is the normalization. 

The prior distribution is found from Jeffreys' rule [17] 

M(f) cc \\P(s \f)[dlnP(s |/)/e/] 2 d5| 1/2 . (11) 

By evaluating numerically the prior distribution (11), and substituting Eq. (8) into Eq. 
(11), and approximating the result by the polynomial of sixth order in / We obtain 

ju{f) = 1.975 - 10.07/+ 48.96/ 2 - 135.6/ 3 

+205.6/ 4 - 158.6/ 5 + 48.63/. (12) 

The distribution P{s\f) takes very small values even for only moderately large 
values of N. Because of this fact, the accurate calculation of the posterior distribution 
becomes a formidable task. In order to simplify the calculation, Eq. (8) may be 
rewritten in the form 

P{s\f) = Q ^(-N </>{/)), (13) 

where 



<f>(f) = (l-f){s) + ^f(0J + 03f)(s 2 ) 

-/in 



l-/+|/(0.7 + 0.V> 



(14) 



Here, the notation 



has been used. One finds the function </> if) to have a deep minimum, say at f=fo- One 
can therefore represent the numerical results in analytical form by parameterizing <f> as 

W) = A +B(f -f ) 2 +C(f -f f, (16) 

where the parameters A,B,C and fo are implicitly defined by (14). 
We then obtain 

P(f \ s )=cf4f)vy(-N[BV -f f+C(f -/ ) 3 ]), (17) 

where c = exp(-NA) / M (s) is the new normalization constant. When P (f\s) is not 
Gaussian, the best-fit value of / cannot be taken as the most probable value. Rather we 
take the best-fit value to be the mean value / and measure the error by the standard 
deviation a of the posterior distribution, i.e. 

f =\[fP(f\sW, and <T 2 =H(f-f) 2 P(f\s)df. (18) 
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The chaoticity parameter / and the standard deviation a are calculated for each group 
of nuclei. The results are given in the following section. 

4. Results 

The search for a phenomenological "control parameter" to describe the evolution 
of the stochastic nature of nuclear dynamics became a subject of nuclear structure in 
the last two decades. 

In the present work, we examine the use of the quadrupole deformation parameter 
as probe of nuclear structure. In spite of absence of complete theoretical study to 
define a certain fixed values of P2 corresponding to critical points of shape/structural 
transitional regions, we tend to classify nuclei into fixed intervals of P2 to get a 
qualitative study. We recall that the analysis of many short sequences of levels tends 
to overestimate the degree of chaoticity measured by a parameter / We focus our 
attention not on the absolute values of /but on the way / changes with p 2 . 

I follow the same method of analysis given in Ref [12]. By grouping nuclei 
according to their P2 deformation parameter, instead of the ratio R4/2 (the ratio 
between lowest 4 + to 2 + level states). The work done in this paper may be a thankless 
work if there is a simple linear relationship between P2 and i?4/2- But on plotting a 
direct plot between them we get figure 2, which shows almost uncorrected relation. 
We start to classify the available even-even nuclei by /?2 to find what can we get? 
Although the capability of the usage of R4/2 to identify the collective excitation of 
nuclei has been confirmed in many published papers depending on both theoretical 
e.g. [18, 19] and empirical [20, 21] calculations. But the R4/2 didn't introduce an 
acceptable method to differentiate between deformed nuclei whether they are oblate 
or prolate deformed, and hence their degree of chaoticity, while the /?2 parameter do 
so. 



-0.05 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 0.50 

Fig. 2: A plot between R4/2 and quadrupole deformation parameter /?2. 

Figure 3 shows a comparison of the spacing distributions that depends on the 
parameter / and the histograms for nuclei divided into classes according to the 
deformation parameter pV In view of the small number of spacings within each class, 
the agreement seems satisfactory. By plotting the best-fit values obtained by Bayesian 
analysis of the chaoticity parameters against P2, we get figure 4. From this figure we 
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see the apparent chaoticity of nuclei having P2 equal or nearly equal to zero these are 
spherical (magic or semi magic) nuclei which are expected to have shell model 
spectra, thus agree with GOE. 
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Fig. 3: NNS distribution of even-even nuclei classified according to deformation 
parameter /? 2 
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Fig. 4: The chaoticity parameter / and standard deviation a (error bars) against 
theoretical deformation parameter fa. 
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An observed minimum of chaoticity is centered at P2 -0.20. Although the 
statistical errors are not good enough to draw a conclusion, but it may indicate that 
something interesting happening through that interval. Two remarkable minima are 
nearly centered at P2 -0.26 and P2 -0.34, these two intervals have nearly the same 
chaoticity parameter that may refer to well deformed nuclei. 

Finally, there is a minimum at P2 -0.14, in order to understand this regular 
behavior we will not study this interval individually, since we observe also a regular 
behavior through its neighbor interval centered at p 2 -0.095. As these two intervals 
are statistically significant from their neighbors, we will extend our study to cover the 
two intervals from p 2 -0.06 to P2 -0.16. By a detailed study of nuclei contributed to 
that interval, we notice nearly an equal number of prolate-oblate nuclei. So, which 
kind is the responsible for this regular behavior? We neglected the negative sign of 
deformation of oblate nuclei, but now the role to test the correctness of such proposal 
comes to surface. We obtained 110 levels of spin-parity 2 + belonging to 15 nuclei of 
oblate shape 62 Ni, 124 Te, 140 Sm, 190 Pt, 192 Pt, 194 Pt, 196 Pt, 198 Pt, 200 Pt, 192 Hg, 196 Hg, 
198 Hg, 200 Hg, 202 Hg, and 204 Hg while 123 levels belonging to 16 nuclei of prolate 
shape 82 Kr, 84 Kr, 88 Kr, 76 Ge, 80 Se, 82 Se, 92 Sr, 96 Mo, 106 Cd, 1?8 Cd, 110 Cd, 112 Cd, 132 Ba, 
136 Ce, 144 Ce, and 192 0s. The NNS distribution is given in figure 5; it is interesting to 
find that the oblate deformed nuclei have more regular spectra than prolate one. This 
may help us to understand the apparent regularity. Indeed, we show [22] that this is 
not only a special case to that interval, but it is the general rule, that oblate nuclei are 
more regular than prolate ones. This fact may be interpreted as, the degree of 
interaction between single particle motion which is chaotic and collective motion of 
whole nucleons which believed to be more regular is weaker in case of oblate 
deformed nuclei than prolate ones. 

Finally, we compare our results based on theoretically calculated P2, with those of 
experimental p 2 values, as deduced from the B(E2; + — > 2 + ) values [10]. Two 
disadvantages of using experimental p 2 arise. One of them is that all nuclei having 
non-zero values, so, it is difficult to differentiate between oblate-prolate deformed 
nuclei. On the other hand, P2 values are not available for all nuclei under investigation 
in this study. We obtained 1050 levels of spin-parity 2 + belonging to 137 nuclei. 
Figure 6 shows the trend of dependence of chaoticity parameter on experimental P2. 
In spite, the error bars are statistical insignificant, the figure shows nearly the same 
trend observed in figure 4. 



5. Conclusion 



We study the nearest neighbor spacing distribution of even-even nuclei classified 
according to the quadrupole deformation parameter using the available experimental 
energy levels. We use a model that interpolates between the Poisson (regular) to 
Wigner (chaotic) distribution by varying the chaoticity parameter from to 1 
respectively. The observed regular behavior of nuclei having certain values of 
deformation parameter indicate that there must be a lot of work to qualify the usage of 
such parameter as an indicator to variation in nuclear structure, and this work may be 
visualized as a step on the road. 
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